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MCUNITV  CLAtUFICATlON 


ADAPTATION  OF  FLUX-CORRECTED  TRANSPORT  ALGORITHMS 
FOR  MODELLING  BLAST  WAVES 


Blast  wave  phenomena  Include  reactive  and  two-phase  flows  associated 
with  the  motion  of  chemical  explosion  products;  the  propagation  of  shocks, 
rarefaction  waves,  and  contact  discontinuities  through  a  nonideal  medium 
(real  air,  possibly  thermally  stratified  and  containing  dust  and  water 
vapor);  and  the  interaction  of  the  blast  waves  (including  boundary  layer 
effects)  with  structural  surfaces.  Flux -Corrected  Transport  (FCT)  represents 
an  accurate  and  flexible  class  of  methods  for  solving  such  nonsteady 
compressible  flow  problems  (Boris  and  Book,  1976).  Coupled  with  a 
nondiffusive  adaptive  gridding  scheme  (Book,  et  al. ,  1980;  Fry,  et  al. , 

1981),  it  enables  complex  time-dependent  shocks  to  be  efficiently 
"captured.” 

In  models  which  treat  all  the  physical  effects  required  for  blast  wave 
simulation,  truncation  errors  inherent  in  the  underlying  finite-difference 
scheme  are  exacerbated  by  nonlinear  coupling  between  the  fluid  equations  and 
by  the  greater  complexity  of  the  phenomena  being  sinulated.  Typical  of  these 
errors  are  the  "terraces"  which  develop  under  some  ci reruns tances  on  the 
flanks  of  sloping  profiles  when  the  growth  of  ripples  due  to  phase  errors  at 
short  length  scales  is  terminated  by  the  action  of  the  flux  limiter.  Two 
approaches  are  possible  toward  eliminating  them:  inproving  the  short- 
wavelength  phase  and  amplitude  properties  of  the  underlying  algorithm,  and 
switching  on  additional  diffusion  locally.  The  latter  approach  folds 
Information  about  the  shape  of  the  profile  and  the  nature  of  the  physical 
process  taking  place  (e.g.,  rarefaction)  Into  the  switch  criterion,  thus 
changing  the  FCT  technique  from  a  "convective  equation  solver"  to  a  "fluid 
system  solver."  In  doing  this,  care  mist  be  taken  to  avoid  losing  the 
accuracy ,  robustness  and  problem-independence  which  constitute  valuable 
attributes  of  FCT  algorithms  (Book,  et  al. ,  198l). 

Tests  carried  out  on  scalar  advection  of  simple  density  profiles  by  a 
uniform  flow  field  show  that  terracing  does  not  require  either  diverging 
velocities  or  discontinuities  in  the  profile,  but  appears  typically  (for  v  > 
0)  where  the  first  and  second  derivatives  of  density  have  the  same  sign 

(Fig.l).  In  order  to  improve  the  properties  of  the  basic  difference  scheme, 

we  propose  a  new  algorithm  for  integrating  generalized  continuity  equations 
over  a  timestep  6t.  Consider  the  following  three-point  transport  scheme: 

Pj  -  p°  -  n(p|+1  -  p°_x)  +  <(p°+1  -  2p®  +  Pj.i); 

Pj  *  pj  -  «ipj+1  -  pj.i)  +  *(pj+i  -  2pJ  +  pj.i); 

pj  =  -  vUj+1/2  "  *J-l/2>» 

where 

♦j+1/2  *  Vl  "  V 


Manuscript  (ubmitted  July  26, 1982. 


The  arrays  {pj }  and  {p ^ }  are  the  old  and  new  densities,  Pj  and  Pj  are 

teiqporary  intermediate  densities,  and  n,  8,  x,  X,  and  p  are  velocity- 
dependent  coefficients*  Here  k  and  X  are  diffusion  coefficients,  and  p  is 


the  antidiffusion  coefficient.  In  the  actual  algorithm,  ♦j+i/g  is  corrected 
(hence  the  name  FCT)  to  a  value  +i/g  chosen  so  no  extrema  in  Pj  can  be 
enhanced  or  new  ones  introduced  in  py  Previous  FCT  algorithms  had  0=0; 

the  widely  used  ETBFCT  and  related  algorithms  (Boris,  1976)  have  in  addition 
<  -  0.  If  we  define  p.  to  be  sinusoidal  with  wave  number  k  on  a  mesh  with 
uniform  spacing  fix,  so'  that  Pj  =  exp  (ij  3)  where  3  =  kfix,  then  the  new 

density  array  satisfies 

pj/pj  s  A  =  l-2i(n+0)sin3  +2(x+X)(cosB-l) 

-  2p(cos3-l)  [l-2insinB+2x(cosB-l) )  * 


From  A  we  can  determine  the  simplification  a  =  A  and  relative  phase  error  R  = 
(l/e3)tan_1(-ImA/ReA)-l,  where  e  =  v6t/6x  is  the  Courant  number. 

Expanding  in  powers  of  $  ue  find  . 

o=l  +  <Xg32  +  ojjB1*  +  og3°  +  .  .  .  ; 

R  =  Rq  +  RgB2  +  +  RgB^  +  .  .  .  . 

First-order  accuracy  entails  making  Rq  vanish,  which  requires  that  n  +  6  = 
e/2.  Second-order  accuracy  (og  =  0)  Implies  that  p  =  x  +  X  -  e2/ 2. 
Analogously,  the  "reduced-phase-error"  property  Rg  =  0  (Boris  and  Book, 

19T6)  determines  p  =  (l-e2)/6,  thus  leaving  two  free  parameters.  One  of 
these  can  be  used  to  make  R^  vanish  also.  The  resulting  phase  error  R(B) 
is  small  not  only  as  3  ♦  0,  but  also  for  larger  values  of  3,  corresponding  to 
the  short  wavelengths  responsible  for  terraces  (Fig.  2).  The  remaining 
parameter  n  can  be  chosen  to  relax  the  Courant  number  restriction  needed  to 
ensure  positivity  from  e  <  1/2  to  e  <  1.  When  coded,  these  changes 
necessitate  a  small  Increase  in  the  operation  count  of  ETBFCT  along  with  a 
Bmall  Increase  in  overhead  to  precalculate  the  two  new  arrays  of  velocity- 
dependent  transport  coefficients.  On  advection  tests,  the  new  algorithm 
completely  eliminated  terraces  (Fig.  3).  When  applied  to  the  coupled  systems 
of  gas  dynamic  equations,  it  produced  profiles  which  closely  approximate  the 
Riemann  solution  of  the  exploding  diaphragm  problem  (Fig.  4). 

The  second  approach  uses  a  rarefaction  flux  limiter  (RFL)  to  eliminate 
numerical  ripples  in  strong  rarefaction  waves.  This  approach  is  physically 
motivated.  Raw  anti-diffusive  fluxes  tj+i/o  are  limited  so  that 
the  slope  of  local  flow  field  profiles  decays  with  time  in  a  rarefaction 
wave.  In  effect,  additional  diffusion  is  left  in  the  field  to  maintain 
monotonicity  of  local  slopes.  For  suit i -material  calculations  a  "contact 
surface  sensor"  is  needed  to  detect  physical  discontinuities  and  shut  off  the 
RFL  locally. 

In  addition  we  found  that  some  care  was  required  when  applying 
generalized  continuity  equation  solvers  to  a  system  of  equations.  Truncation 
errors  of  the  various  equations  can  interact,  causing  undershoots  or 
overshoots  in  nonconvective  quantities  such  as  pressure.  We  found  that  it 
was  necessary  to  monotonlze  derived  quantities  (pressure,  velocity)  before 
using  them  in  minimal-diffusion  transport  algorithms. 


The  above  methodology  has  been  applied  to  a  series  of  test  problems 
initiated  tor  a  spherical  high-explosive  (HE)  detonation  in  air.  An  ideal 
Chapman -Jouguet  detonation  was  used  to  specify  the  initial  conditions; 
afterburning  was  neglected.  In  the  absence  of  reflecting  surfaces, 
spherical  symmetry  is  maintained  and  the  calculation  remains  one- 
dimensional.  A  nonuniform  radial  grid  vas  used  with  extremely  fine  zoning 
near  the  shock  front.  The  grid  vas  moved  so  that  the  shock  remained 
approximately  fixed  with  respect  to  the  mesh.  The  original  version  of  the 
FCT  algorithm  gave  rise  to  pronounced  terraces  in  the  rarefaction  region. 

This  would  have  rendered  any  two-dimensional  calculation  involving  shock 
diffraction  or  nonideal  effects  dubious.  The  techniques  described  here 
inproved  the  blast  wave  results  considerably.  The  decrease  in  phase  error 
reduced  terracing  dramatically. 

Next,  a  two-dimensional  (2D)  numerical  calculation  vas  performed  to 
simil&te  one  of  Carpenter's  ( 197*0  height -of-burst  experiments  which  used 
spherical  8-lb.  charges  of  PCX  9**0U  at  51*6  cm.  The  previous  fine-zoned  ID 
calculation  was  used  to  initialize  the  problem.  It  was  napped  onto  the  2D 
grid  Just  prior  to  the  onset  of  reflection.  The  solution  was  then  advanced 
in  time,  with  pressure  being  calculated  from  a  real-air  equation  of  state  and 
a  JWL  equation  of  state  for  the  combustion  products.  The  front  of  the  blast 
wave  was  captured  in  a  finely  gridded  region  vhich  moved  outward  horizon¬ 
tally.  Special  care  vas  taken  to  ensure  that  the  grid  moved  smoothly.  The 
resulting  solution,  particularly  the  curve  of  peak  overpressure  vs.  range, 
was  consistent  with  Carpenter's  experimental  data  (Fig.  8).  Although  this 
calculation  represents  a  reasonable  accurate  sinulation  of  the  double-Mach- 
8 tern  region,  no  doubt  improvements  can  and  will  be  made  to  numerically  model 
such  phenomena. 
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Fig.  1  —  Rounded  half  circle  used  in  passive  scalar  advection  tests  (a)  initially, 
and  (b)  after  propagation  for  14  cycles  using  JPBFCT.  Note  that  terraces  form 
even,  as  here,  in  the  absence  of  comers  in  the  profile.  Tick  marks  indicate 
computational  zones  (N  *  100). 


Fig.  2  —  Contour  plot  of  R(0,e)  for  new  multicoefficient  FCT  algorithm. 
Note  R  %  0  except  for  6  >  3  it  12.  The  relative  phase  error  vanishes  ex¬ 
actly  for  e  =  1/2  and  e  ■  1 . 


JPBFCT  ROUNDED  HALF-CIRCLE 
CYCLE  -  1*  TRUE  -  2.17  10  * 


LCPFCT  ROUNDED  HALF-ORCLE 
CYCLE  -14  TIME -2.17  10 


(b) 

Fig.  3  —  (a)  Blowup  of  Fig.  1(a)  (dashed  line)  compared  with  (b)  same  profile  as 
computed  using  new  sixth  order-phase  -accurate  FCT  algorithm.  Solid  traces  are 
exact  solutions. 


HE 
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-AIR  JPBFCT  S  T «  1.00*10 
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PRE55URE  VS  TIME 

Fig.  7 — Pressure-time  histories  directly  beneath  burst  site.  Note  second  peak,  associated  with 
interaction  between  shock  reflected  from  ground  and  following  contact  surface. 
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01CY  ATTN  DRCDE’O  L  FLYNN 

01CY  ATTN  DR  YAM- Tl  (TECH  LIB)  UNCL  ONLY 


CCMMANDEP 

U  S  ARVY  MISSILE  C.CMNA*0 
PEOSTONE  ARSENAL,  AL  9 

OICY  ATTN  PSIC 
01CY  ATTN  OR.OVI-XS 


COMMANDER 

II  S  ARMY  MOBILITY  SCLIR  PCD  CMO 
FORT  BELVOIR,  VA  22060 

(CNWDI  TO  APMY  MAT  OEV  £  PEADINESSS  COMMANDl 
OICY  ATTN  DROME-VC  (TECH  LIB) 

COMMANDER 

U  S  ARMY  WUCLF4P  £  CFFMICAL  AGENCY 
7500  9ACKLICK  RCAQ 
B?JILD TNG  2073 
SPRINGFIELD,  VA  221 50 

(DESIRES  ONLY  1  CY  TO  LIBRARY) 

OICY  ATTN  J  SIMMS 
OICY  ATTN  LIPPARY 


COMMANDANT 
IJ  S  ARMY  WAR  COLLEGE 
CAPLISLE  BARRACKS,  PA  17013 
OICY  ATTN  LIBRARY 

C?VNAMDFP 

DAVID  TAYLD?  NAVAL  SUP  R  f.  D  CTP 
EETHESDA,  mq 

ICNWOT  CNIY  ATTN  ».  EIPKHFAD  CCDE  5915.6) 

OICY  ATTN  CODF  142-3  (UePAP.Y) 

CFFICER-IN-CMARGE 

NAVAL  CIVIL  ENGINES PING  LA9CPAT0RY 
PORT  HHFNPMS,  CA  030AI 

OICY  ATTN  COOE  L5?  J.FCRPFST 
OICY  ATTN  CODF  LCPA  (LIBRARY) 

OICY  ATTN  CODE  L51  J  CRAWFORD 

OICY  ATTN  L 51  R  MLRTHA 

# 

COMMANDER 

NAVAL  ELECTRONIC  SYS7FMS  COMMAND 
WASHINGTON,  DC  2036C 
OICY  ATTN  OME  117-21 
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COMMANDER 

NAVAL  FACILITIES  ENGINEERING  COMMAND 
WASHINGTON*  DC  203S0 
.  OICY  ATTN  CODE  04e 

HEADQUARTERS 
NAVAL  MATERIAL  COMMAND 
WASHlNGTON«  DC  20‘*60 
OICY  ATTN  MAT  08T-2? 

commander 

\AVAl  OCEAN  SYSTFW$  CFMT£° 

SAN  DIEGO*  C A  *?1S? 

DlC Y  ATTN  CODE  013  E  CCCPFO 
OICY  ATTN  CODE  4471  I  TECH  LID) 


SUPERINTENDENT 

NAVAt  POSTGRADUATE  SCHOOL 

MCNTFREV,  CA  03940 

f DESIRES  NC  CNKCI  CCCUMENTS) 

OICY  ATTN  CODE  1424  LIBRARY 
OICY  ATTK  C-  LINDSAY 

COMMANDING  OFFICcR 
NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  DC  2C37S 

CRD  r.  pD/m  ATTN  COPE  1221  F!?R  £  PPD  ATTN  CODE  2628  FOP) 
OICY  ATTN  CODE  40*0  J  *0*1? 

20CY  ATTN  CPDF  2627  I  TECH  LIB) 

100CY  ATTN  CODE  4C4C  C  3rrK 

COMMANDER 

NAVAL  SEA  SYSTEMS  CCMMAKO 
WASHINGTON,  DC  20362 

OICY  ATTN  SEA-09GS?  CLIP) 

OICY  ATTN  SRA-03 SI 

DPPJCER  IN  CHARGE 
NAVAL  SURFACE  WfARCKS  C5NTFP 
WHI TE  CAK  l AR ORATORY 
SUV=R  SPRING*  *D  2CSK 
OICY  ATTM  P44  H  CLAZ 
DIC Y  ATTN  CODE  F31 
OICY  ATTN  CrOE  Y211  (TECH  LIB) 


xovmanoer 

NAVAL  SUPEACF  WFAR.7KS  CFNTFP 
DAHLGPEK*  VA  22*4« 

D1CV  ATTN  TECH  LIBRARY  £  INFO  SVCS  SR 


IS 


PRESIDENT 
NAVAL  MAR  CCLIEGF 
NEWPORT,  »I  0 7R40 

01CY  ATTN  CODE  E-ll  (teCH  SERVICE) 
COMMANOE® 

naval' WEAPONS  CENTER 
CHINA  LAKE,  CA  R3S5* 

OICY  ATTN  CODE  3  201  P  CrpCLE 
01CY  ATTN  C^DE  266  C  AUSTIN 
OICY  ATTN  COOP  22?  (TECH  L!0) 

COMMANDING  OFFICER 
NAVAL  WEAPONS  EVALUATION  FACILITY 
KIRTLANO  AIR  FORCE  EASE 
ALBUOUERC'JE ,  NM  P7117 
OICY  ATTN  R  HUGHES 
OICY  ATTN  CODE  1C  (TECH  LIB) 

OFFICE  OF  NAVAL  FESEAPCF 
ARLINGTON,  VA  22217 

OICY  ATTN  CODE  474  N  PERPPNP 


•r ettire  rr  tjjc  rgjcp  pc  NA^AL  PRFRATICVS 

KASHf  \'G-rr»  ,  20?5C 

DIC Y  ATTN  SB  L 
OICY  ATTN  C3  C3EG 

DIRECTOR 

STP  AT  EGTC  SYSTCMS  PPCJECT  CFF Ire 
DE3APTMgNT  OF  THT  *«AVY 
WASHINGTON,  DC  20376 
OICY  ATTN  NS P-272 
OICY  ATTN  NSP-43  (TECH  LIP) 

.•AIR  PC’JCE  GE^*»HVSICS  L  ABCRA  Troy 
HANSCOM  A rBf  vA  a1711 

OICY  ATTN  LWW  K  THrvcSON 

AIR  FORCE  INSTITLTE  OF  TECHNOLOGY 
AIR  UNIVERSITY 

WRIGHT- PATE* SDN  AFP,  OH  4543? 

(DOES  NOT  DESIRE  CLASSIFIED  DOCUMENTS) 
OICY  ATTN  LIBRARY 

HCADOUAPTEPS 

AIR  FORCE  SYSTEMS  CCHMANC 
AMOREWS  AFB,  OC  2033 A 
OICY  ATTN  DLWM 
OICY  ATTN  DLW 


AIR  FORCE  WEAPONS  LAPCPATnpv,  A  ESC 
KIPTLANO  AF9,  N*  87117 


OICY 

ATTN 

NTFS-C  R  HENNY 

OICY 

ATTN 

NTEO-I 

OICY 

ATTN 

NTED  R  KATALUCC! 

OICY 

ATTN 

NTF  "  PLAMCAOCN 

OtCY 

ATTN 

NT  D  PAYTCN 

OICY 

ATTN 

NTED- A 

OICY 

ATTN 

NTES-G 

OICY 

ATTN 

SUL 

OICY 

ATTN 

OEY 

OICY 

ATtnj 

NTFS-S 

OICY 

ATTN 

NTFC 

OICY 

ATTN 

DEV 

DI°ECT0° 

AIR  UNIVEPSITv  LIBRARY 

ncoAPTwtNT  OF  THF  AIP  FCPCE 
MAXELL  AFR ,  AL  16112 

(HESJPFS  KC  CNWD1) 

OICY  ATTN  AUL-LSE 

ASSISTANT  CHIEF  OF  STAFF 
INTELLIGENCE 

DEPAPTMEMT  r»F  THE  AIR  FORCE 
WASHINGTON.  DC  2033C 

OICY  ATTN  IN  PM  4AS32 

* 

ASSISTANT  CHIEF  CF  STAFF 
STUDIES  £  ANALYSES 
PE°A° TMENT  OF  THE  AIR  FORCE 
WASHINGTON,  OC  20320 

OICY  ATTN  AF/5.A*I  ( TFCH  119) 

ASSISTANT  SECRFTAOY  OF  THE  AF 
BESEARCH,  DEVELOPMENT  C  LOGISTICS 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  OC  2033C 

OICY  ATTN  SAFALP/CEP  FOR  STRAT  £  SPACE  SYS 


BALLISTIC  MISSILF  1FFICE/MN 
A I o  FOPCE  SYSTEMS  CC"EA*T 
NORTC*'  AFB,  CA  R?ACR 
triMUTCNANI 

OICY  ATTN  »NNXH  G  KALANSKV 

OICY  ATTN  MNNXH  W  CFLVECCHIO 

OICY  ATTN  MNN  W  CPADTPEE 

OICY  ATTN  MNMXH  D  GAGE 

OICY  ATTN  MMMX 


SO 


'1 

•.1 

,  <i 

N 


r 


OE°UTY  CHIEF  OF  STA'F 
RESEARCH,  DEVELOPMENT,  C  ACC 
OEPARTMENT  OF  TWF  A  IF  FCFCF 
WASHI NGTCM «  DC  ?0?3Q 

OICY  ATTN  AFFCOI  K  ALFXANO»OW 
OICY  ATTN  A«*RPON 
OICY  ATTN  AFPDGl 

OE»U TY  CHIEF  OF  STAFF 
LOGISTICS  6  ENGINEERING 
DEPARTMENT  OF  THE  A  I®  FORCE 
WASHINGTON  »  OC  2033C 
OICY  ATTN  L6FF 


COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION*  AFSC 
WRIGHT-oaTTEPSON  AFE,  CH  45433 
OICY  ATTN  Ml  IS  L IPRARY 


COMMANDER 

povp  AIR  DEVELOPMENT  CENTER,  AFSC 
GPTFFISS  AF3.  NY  13441 

ICESI"FS  NT  CNWDT) 

OICY  ATTN  TSLO 


B 


r  . 

r- 

r, 


* 

'1 


STR  ATEGIC  A  IP  COMMAND 
DEPARTMENT  OF  THF  AI»  F&PCE 
OFFUTT  AFB,  NB  6M13 

OICY  ATTN  NPI-STINFC  LIBRARY 
OICY  ATTN  XPFS 

OICY  ATTN  I NT  J  MCKINNEY 

* 

VELA  SEISMOLOGICAL  CENTEP 
312  yONTGrwER Y  STREET 
ALEXANOIRA,  VA  22314 
OICY  ATTN  G  '.'ILRICH 


OE°AP TMf  NT  OF  ENERGY 
ALSUCUEPOijc  PPCPATICNS  CFFICF 
P  0  BOX  5400 
AL9U0'JF°0UE  *  NN  P711B 
OICY  ATTN  CTID 


R  OEPART«FNT  op  ENFPGY  ■«- 

-J  WASHINGTON,  oc  2054  5 

OICY  ATTN  CMA/RC?67 

DEPARTMENT  OF  ENEPGY 
*•  NEVADA  OPERATIONS  OFFICE 

M  P  0  SOX  14100 

7  LAS  VEGAS,  NV  99114 

OICY  ATTN  MAIL  £  PECOROS  FOR  TECHNICAL  IJRPARY 

V 

i 
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LAWPENCF  LIVERMORE  NATIONAL  LAB 
P  0  box  eop 
LIVEPMOREf  CA  04550 
OICY  ATTN  L**°0  R  OCMG 

OICY  ATTN  1.-205  J  FEAPST  (CLASS  L-203) 

OICY  ATTN  L-90  9  NOPRIS  (CLASS  L-504) 

OICY  ATTN  L ~ 7  J  KAHN 

OICY  ATTN  0  GLENN 

OICY  ATTN  L  437  P  SCHCCK 

OICY  ATTN  TECHNICAL  INFO  D*?T.  LIBRARY 

OICY  ATTN  L-200  T  BLTK0V1CH 

LOS  ALAMOS  NATIONAL  SCIENTIFIC  LAB 
''AIL  ST  AT  ION  500  0 
P  n  BOX  1463 
LOS  ALAMOS «  MM  P7545 

(CLASSIFIED  ONLY  TC  MAIL  STATION  50001 
OICY  ATTN  R  WFITTAKFR 
OICY  ATTN  C  KELLER 
OICY  ATTN  MT  SANFORD 
OICY  ATTN  MS  364  (CLASS  REPORTS  LIBI 
OICY  ATTN  E  JONES 


LOVELACE  BIOMFOICAL  € 

ENVIRONMENTAL  PSCH  INSTITUTE,  INC* 

P  0  BOX  5690 
ALSUOIJFPgiJF,  NM  P7115 

OICY  ATTN  °  JCN'S  OJNCI  CN’LY) 

OAK  PIHGE  NAT  10*! AL  LABORATORY 
NUCLEAR  DIVISION 
X-10  LAB  RECORDS  DIVISION 
P  0  BOX  X 

CAK  PIDGE,  TN  37830 

OICY  ATTN  CIVIL  DEF  RFS  PPOJ 
OICY  ATTN  CENTRAL  PSCH  LIBRARY 

SAMOIA  LABORATORIES 
LIVERMORE  LABORATORY 
o  0  POX  969 
LI  vcp.MrjpE  ♦  CA  94559 

OICY  ATTN  LIBRARY  C  ScCU«>ITv  CLASSI  r  IC  AT  I  CM  DIV. 


SA*‘DIA  NATIONAL  LAB 

p  o  pox  5eoo 

AI.BUOUEROUS,  N‘»  S71  95 

(ALL  CLASS  ATTN  SCC  CDNTPQL  OFC  FOR) 
OICY  ATTN  A  CHA3AN 
OICY  ATtn  t  HILL 
OICY  ATTN  0?C  1  ? 50  W  BROWN 
OICY  ATTN  A  CI-ABIA 
OICY  ATTN  W  PCHERTY 
OICY  ATTN  314! 

OICY  ATTN  L  VGPTNAM 
OICY  ATTN  J  BENISTER 

CENTRAL  I  NT  ELL  I GFMCE  AOEFCY 
WASHINGTON,  OC  205C5 
OICY  ATTN  OSKP/NEO 

DEPARTMENT  OF  THF  INTER ICP 
BUREAU  CF  MINES 

BLDG  20,  OENVE9  FEDERAL  CENTFP 
DENVER,  CO  P0225 

I (UNCI  CNL Y )  | 

OICY  ATTN  TECH  LIP  (UKCl  CN'LYl 
OtPEC TOR 

FEDERAL  EMERGENCY  MANAGEMENT  AGENCY 
NATIONAL  SEC  OFC  LITIGATION  f.  RSCH 
1725  I  STREET,  NW 
WASHINGTON,  OC  20472 

((ALL  CLASS  ATTN  BIC5  OCC  CCNTPCL  FCR) ) 
OICY  ATTN  MITIGATION  5  RSCH  DIV 

AC'toFX  C09P. 

4«5  CLYDF  AVFVljr 
MOUNTAIN  VIEW,  CA  <*4042 
OICY  ATTN  C  WfLF 

A EROS  RACE  COPP. 

P  C  BOX  92957 
LOS  ANGELFS,  CA  90009 
OICY  ATTN  H  LIPSLS 

OICY  ATTN  TECHNICAL  INFORMATION  SERVICES 

AGP  AB I AN  ASSOCIATFS 
250  N  NASH  STREET 
El  SFGUNDO,  CA  90245 
OICY  ATTN  »  AGBABIAN 

ANALYTIC  SERVICES,  INC. 

400  ARMY-NAVY  OR TVE 
ARLINGTON,  VA  22202 

OICY  ATTN  G  HESSFLBACHSR 


APPLIED  PESEA°CW  ASSOCIATES,  INC 
?601  WYOMING  PLVH  NE  S''iTE  F-l 
ALBUQUERQUE  ,  N”  P711? 

01CY  ATTN  J  FCATTCN 
0 1  CY  ATTN'  N  HIGGINS 

APPLIED  THEORY.  INC. 
mio  west  upon  ?lvd 
ms  ANGELES *  CA  50^24 

<?  CVS  IE  LNCLASS  PP  l  CY  IF  CLASS ) 
OICY  ATTN  J  TRULIC 

AOTEC  ASSOCIATES,  INC. 

26046  EDEN  LANDING  PPAD 
HAYWARD,  CA  94545 
OICY  ATTN  S  GILL 

4 STPPN  FESAPCH  £  ENGINEERING 
1001  OLD  MI COLFF IOLD  NAY  #15 
MOUNTAIN  VIEW,  CA  54043 
OICY  ATTN  J  HUNTINGTON 

AVCO  RESEARCH  C  SYSTEMS  CROUP 
201  LOWE LI  STREET 
WILMINGTON,  MA  918*7 

OICY  ATTM  LIBRAPY  AF30 
* 

BOV  COPP. 

7915  JONES  BPANC.H  DRIVE 
VC L  FAN,  VA  ??1Q2 

OICY  ATJN  ft  LAVAGNINC 
OICY  ATTN  T  MEIOHEfPS 
OICY  ATTN  CORPORATE  LIERAPV 

BOM  rnpo, 

°  0  BOX  9274 
AL BUO'JEP DUE ,  Nv  87119 
OirY  ATTN  R  HENSLEY 

BOEING  CO. 

P  0  BOX  3707 
SEATTLF,  WA  9R124 
OICY  ATTN  S  $T° AC  K 
OICY  ATTN  AFROS0  ACF  LIBRAPY 
OICY  ATTN  M/S  42/37  P  CA°l$ON 

CALIFORNIA  PESEAPCH  S  TECHNOLOGY,  INC. 

6269  VAPTEL  AVENUE 
WOOD!  AND  HILLS,  CA  91367 
OICY  ATTN  LTBPARV 
OICY  ATTN  K  KREYENHAGEN 
OICY  ATTN  M  ROSENBLATT 


CAL  1 FPRN  T  A  RESEARCH  f.  TECHNCLOGY,  INC. 
4049  FIRST  STRFFT 
L I V CF  MD R F »  CA  c4550 
OICY  ATT N'  3  PRPHAL 

C  Al  SD  AN  COR®. 

®  0  POX  4r’0 
BUFFALO,  NY  14275 
0  ICY  ATTN  UP»AnY 

OENVRR,  UNIVERSITY  OF 
COL  CP  A0P  SEMINipv 
OPNVFR  pFSFAPCH  I N $ T  I T Ut 0 
r>  0  POX  12177 
OE'.'VFR,  CO  80210 

(ONLY  1  COPY  Op  CLASS  RPTS) 
01CY  ATTN  SEC  OPFICFR  FCR  J  WISOTSKI 

EGtG  WASH.  4NAI  YTICAL  SVCS  CTP,  INC. 
o  o  pry  102 1 P 
ALP'JOUEROUE,  NM  67114 
0! Cv  ATTN  LIEF ARY 

FRTC  H.  WANG 

CIVIL  ENGINEERING  DSCH  FAC 
UNIVERSITY  OP  NEK  MEXICO 
UNIVFRSITY  STATirv 
P  0  20X  25 

ALBUGHERQME  •  NY  87121 
0 1C Y  ATTN  J  LAMB 
OICY  ATTN  P  LCO0E 
OICY  ATTN  N  BAUM 
OICY  ATTN  J  KOVAPNA 

GA°0»  INC. 

744°  N  NATCHEZ  AVCNLE 
Nil c S  t  TL  60646 

01CV  ATTN  G  NEIGHArOT  fif'd  ONLY ) 

GENE®  AL  ELECTRIC  CO. 

SPACE  DIVISION 
VALLEY  FPRGE  SPACE  CENTER 
D  0  enx  9555 
PHI LAOELPHI A»  PA  19101 
OICY  ATTN  M  30PTNFP 

GENERAL  RESEARCH  COPP. 

SANTA  PARRARA  DIVISICN 
P  0  BOX  6770 
SANTA  BARBARA,  CA  92111 
OICY  ATTN  TTO 
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M-TFCH  labs*  INC. 

P  0  BOX  1686 
SANTA  MONICA*  CA  <>0406 
01CY  ATTN  3  HAPTENC ALP 

HORIZONS  TECHNOLOGY,  INC. 

7330  CLA1PE«0NT  KESA  BLVO 
SAN  0 1 E GO ,  CA  92111 
01 CY  ATTN  R  KPllGEP 

I I T  RESEARCH  INSTITITE 
10  W  ?5TH  STREET 
CHICAGO.  TL  60fe ! 6 
0 1C Y  ATTN  P.  WELCH 
OICY  ATTN  M  JOHNSON 
01CY  ATTN  DOCUMENTS  LIBRARY 


INFC°MAT ION  SCIENCE,  INC. 
123  W  f*A9RF  ST°£FT 
SANTA  BARBARA  ,  CA  93105 
OICY  ATTN  W  01107  I AK 


INSTITUTE  FOP  DEFENSE  ANALYSES 
400  APMy-MAVY  Dn IVE 
ARLINGTON,  VA  22202 

OICY  ATTN  CLASSIFIED  LIBRARY 

J  D  HALTIWANGER  CCNSLLT  EN6  SVCS 
RM  106A  CIVIL  ENGINEERING  BLOC 
208  N  ROMINF  STPE6T 
••RRANA,  IL  61*01 
OICY  ATTN  W  HALL 

J.  H.  WIGGINS  CO.,  INC. 

1650  S  PACIFIC  COAST  HIGHWAY 
REDONOH  BEACH,  CA  SC277 
OICY  ATTN  J  COLLINS 

KAMAN  AVI9YNF 
83  SECOND  AVENUE 
NO»THWEST  IN0USTP1AI  PARK 
B’J®  LI  MGTCN  ,  MA  01*03 
OICY  ATTN  P  PLETENIK 
OICY  ATTM  LIBRARY 
OICY  ATTN  N  HOBBS 
OICY  ATTN  E  CRISCICNE 

KAMAN  SCJFNCFS  CO P* . 

R  0  BOX  7463 

COLORADO  SPRINGS,  CC  80933 
OICY  ATTM  D  SACHS 
OICY  ATTM  F  SHELTCN 
OICY  ATTN  LIBRARY 
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KAMA*  TEMPO 

916  STATE  STREET  <o  0  DRAWER  00) 

SANTA  BARBARA,  CA  92102 
0 ICY  ATTN  DASIAC 

LOCKHEED  MISSILES  6  SPACE  CC.  ,  INC. 

P  0  POX  504 
SUNNYVALE,  CA  94086 
01CY  ATTN  J  WEISNER 
OICY  ATTN  Tlf-liePARY 

MARTIN  MARIETTA  COR®. 

P  0  BOX  5827 
OP l ANDO,  EL  32855 
OICY  ATTN  G  FCTIFC 

MARTIN  VAP I  ETTA  fORp. 
o  0  BOX  179 
DENVER,  CO  R0?01 
OICY  ATTN  G  FREVER 

MCDCMNELL  CC'.'GLAS  CCPP. 

5131  PT>LSA  AVENGE 
HUNTINGTON  95 ACH,  CA  92647 
OICY  ATTN  H  ME® OMAN 
OICY  AT^N  R  HALPRIN 
OICY  ATTN  0  OF AN 

MCOONNELL  DOUGLAS  CO®P.  , 

1955  LAKEWOOD  BOULEVARD 
LONG  BEACH,  CA  90846 

OICY  ATTN  M  POTTER 

• 

MERRITT  CASES,  INC. 

P  0  BOX  1206 
REDLANDS*  CA  92373 
OICY  ATTN  J  MTRRITT 
OICY  ATTN  LIPPARY 

METEOROLOGY  RESEARCH,  INC. 

464  W  WOCnBU°Y  RCAO 
ALTADEMA,  CA  91001 
OICY  ATTN  W  GREEN 

MISSION  PESEARCH  CORP. 

P  0  OPAWF*  719 

SANTA  BARBARA,  CA  93107 

(ALL  CLASS:  ATTN:  SEC  OFC  FOP) 
OICY  ATTN  C  LONGMIPE 
OICY  ATTN  G  MCCARTCR 
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**4C  IF  TC- SI  ESP.  A  RESEARCH  CHPP. 
’1456  CL  OVER  FI  flf?  BLVC 
SANTA  MONICA,  CA  904C4 
OICY  ATTN  M  3P00E 


PACIFIC-SIERRA  PESEARCH  CHPP. 

WASHINGTON  OPERATIONS 
1401  WILSON  SU'D 
SUITE  1100 

ARLINGTON,  VA  22209 
01CY  ATTN  0  GORMLEY 

PACIFICA  TECHNOLOGY 
P  0  BPX  149 
DEL  MAR,  CA  92014 
01C  Y  ATTN  R  P  JCRK 
01CY  ATTN  G  KENT 
OICY  ATTN  TETH  LIPRAPY 

PATFL  ENTERPRISES,  INC. 

P  0  BOX  3531 
HUNTSVILLE,  AL  3591C 
OICY  ATTN  «  PATEL 

PHYSICS  INTERNATIONAL  CO. 

2700  MERCED  STREET 
SAN  LEANOPO,  CA  94577 
OICY  ATTN  L  PEHRMANN 
OICY  ATTN  TECHNICAL  LIBDAPY 
OICY  ATTN  E  MOORE 
OICY  ATTN  j  T HCV S£*: 

OICY  ATTN  F  ^AIJFP 

R  C  0  ASSr'CIATrs 
p  0  POX  °655 

MAP IN A  CEL  RET,  CA  9C291 
OICY  ATTN  R  RHPT 
0 ICY  ATTN  A  KUHl 
OICY  ATTN  J  L  P'/.  !  S 
OICY  ATTN  W  WRIGHT 
OICY  ATTN  J  C A° ° E NT E P 

OICY  ATTN  TECHNICAL  INFORMATION  CENTER 


PAND  COPP. 

1700  ^AIN  STREET 
SANTA  *CNICA,  CA  ^0406 
OICY  ATTN  C  *CW 

SCIENCE  APPLICATIONS,  INC 
RADIATION  INSTRUMENTATION  OlV 
4615  HAWKINS,  Nf 
ALBUO'JFRO'JE,  NY  9 Tl CO 
OICY  ATTN  J  OfSHCN 
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SCIENCE  APPLICATIONS,  INC. 
n  0  OOX  2351 
14  jnilA,  C 4  9203? 

OICY  ATTN  w  X'lLSCN 

01CY  ATTN  TECHNICAL  LIBRARY 

OICY  ATTN  R  SCHLAU C 

SCIENCE  APPLICATIONS,  INC. 

101  CONTINENTAL  8LV0 
FL  SEG»»NOO,  CA  OQ24  5 
OICY  ATTN  0  HCVE 

SCIENCE  APPLICATIONS,  INC. 

2450  WASHINGTON  AVENUE 
SAN  LEANDRO,  CA  94577 
OICY  ATTN  0  BERNSTEIN 
OICY  ATTN  D  MAXWELL 

SCIENCE  APPLICATIONS,  INC. 

P  0  POX  1303 
*CLEAN,  VA  ??l 02 

OICY  ATTN  J  COCKAYNE 
OICY  ATTN  B  CWAMBEPS  III 
OICY  A^TN  M  K NAS El 
OICY  ATTN  W  LAYSCN 
OICY  ATTN  R  SIEVERS 

SOUTHWEST  RESEAPCH  INSTITUTE 
»  C  ORAWER  78510 
SAN  ANTCNIfl ,  TX  78284 
OICY  ATTN  A  WENZEL 
OICY  ATTN  W  BAKER 

SRI  lNTEPNATICNAL 
333  PAVENSKOrD  AVFNUE 
MENLO  PARK,  CA  94075 

OICY  ATTN  G  APRAHAVSCK 
OICY  ATTN  LIPPARY 
OICY  ATTN  J  CCLTCN 

SYSTFMS,  SCIENCE  S  SOFTWARE  INC 
P  0  BOX  8243 
AIBUOUERQHF  wm  P7t08 
OICY  ATTN  C  NEEDHAM 


SYSTEMS,  SCIENCE  C  SOFTWARE,  INC. 
P  0  BOX  16?0 
LA  JOLLA,  CA  92038 
01CY  ATTN  J  BARTWEL 
0 1C Y  ATTN  T  R  INF Y 
OICY  ATTN  0  GRINE 
OtCY  ATTN  LI  BPARY 
OICY  ATTN  C  HASTING 
OICY  ATTN  K  PVATT 
OICY  ATTN  C  OISMliKFS 
OICY  ATTN  T  CHERRY 

SYSTEMS,  SCIENCE  £  SOFTWARE,  INC. 
llo00  SUNRTSF  VALLEY  CRIVE 
RESTCN,  VA  22^91 
OICY  ATTN  J  MURPHY 


TSLEOYNE  OPOWN  ENGINEERING 
CimmtnGS  PESEAOCH  »APK 
HUNTSVILLE,  AL  3**07 

OICY  ATTN  J  PAVENSCRAFT 
OICY  ATTN  J  MCRWAJN 


TE*°A  TEK,  INC. 

420  WAKAP.A  WAY 
SALT  LAKE  CITY,  LT  E41C9 
OICY  ATTN  A  ABOM-SAVEC 
OICY  ATTN  Lie«>AP.Y 
OICY  ATTN  A  JONES 
OICY  ATTN  S  GREEN 

TETRA  TECH,  INC. 

830  N  ROSEMEAD  BLVO 
PASADENA,  CA  91107 
01C\  ATTN  L  HWANG 

TRW  DEFENSE  £  SRACF  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  «0?79 
OICY  ATTN  N  LIPMER 

OICY  ATTN  TECHNICAL  INFORMATION  CENTER 
OICY  ATTN  T  MAZZCLA 

TP.W  DEFENSE  £  SPACE  SYS  GROUP 
o  O  BOX  1310 

SAN  BERNARDINO,  CA  S240? 

OICY  ATTN  G  HL'LCHEP 
OICY  ATTN  P  DAI 
OICY  ATTN  E  WCNG 


UNIVFPSAL  ANALYTICS,  I*C. 
7740  W  MANCHESTER  RLVC 
**L AYA  OEt  REV,  C*  0C?«?i 
OICV  ATTN  E  F IELO 


WEI 01  INGEP  ASSOC.,  CCNSULTING  ENGINEERS 
IIO  E  5<»TH  STREET 
NEW  YORK,  MY  10072 
OICV  ATTN  I  SANOLEP 
OICY  ATJN  M  PAPON 


NEIDIINGER  ASSOC.,  CCNSLLTING  ENGINEERS 
3000  SAND  HILL  PCAD 
MENLO  PARK,  CA  04075 
OICY  ATTN  J  I SSNPEPG 


